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Abstract 

Dipole models based on various saturation scenarios provide reasonable fits to small-x DIS 
inclusive, diffractive and exclusive data from HERA. Proton un-integrated gluon distributions 
extracted from such fits are employed in a fc_L -factorization framework to calculate inclusive 
gluon distributions at various energies. The n-particle multiplicity distribution predicted in the 
Glasma flux tube approach shows good agreement with data over a wide range of energies. 
Hadron inclusive transverse momentum distributions expressed in terms of the saturation scale 
demonstrate universal behavior over a wider kinematic range systematically with increasing 
center of mass energies. 
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1. Introduction 

HERA deeply inelastic scattering (DIS) results on structure functions demonstrate 
a rapid bremsstrahlung growth of the gluon density at small x. When interpreted in 
the same framework as the parton model, this growth is predicted to saturate because 
the gluon occupation number in hadron wave functions saturate at a value maximally 
of order 1/as', dynamically, nonlinear effects such as gluon recombination and screen- 
ing by other gluons deplete the growth of the gluon distribution[l]. Gluon modes with 
kx < QsO^ ^qcd) are maximally occupied, where Q'g{x) is a dynamically generated 
semi-hard scale called the saturation scale. For small x, Qs{x) is large enough that high 
occupancy states can be described by weak coupling classical effective theory [2]. This 
Color Glass Condensate description of high energy hadrons and nuclei is universal and 
has been tested in both DIS and hadronic collisions. In particular, saturation based 
phenomenological predictions successfully describe recent LHC p+p data [5,6] and pre- 
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diet possible geometrieal sealing of transverse momentum distribution[4,7] similar to the 
geometrieal scaling observed previously in DIS. 

The object common to DIS and hadronic collisions is the dipole cross section -grg^ ('"J- ' -^i bj^). 
In the CGC framework, the dipole cross section can be expressed in terms of expectation 
values of correlators of Wilson lines representing the color fields of the target. The energy 
dependence of this quantity comes from renormalization group evolution but to get the 
realistic impact parameter dependence one has to rely on models involving parametriza- 
tions constrained by experimental data. In the large Nc limit, the dipole cross section is 
related to the un-integrated gluon distribution inside hadron/nucleus as 
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For hadron-hadron collisions, the inclusive gluon distribution which is k^-factorizable 
into the products of un-integrated gluon distributions in the target and projectile is 
expressed as 
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2. Saturation models of HERA DIS 



Two models of the dipole cross-section that have been extensively compared to HERA 
data are the IP-Sat [8,9] and the b-CGC [10,11] models. In the former the impact pa- 
rameter dependence is introduced through a normalized Gaussian profile function Tp(h±) 
and in the latter through a scale Qs{x,h±). For a detailed discussion of the parameters 
involved in these models and their values from fits to HERA data, see ref. [6]. The sat- 




Fig. 1. Left: Adjoint saturation scales obtained from the IP-Sat(blue lines) and b-CGC models(red 
crosses). Right: Impact parameter and energy dependence of Q^Sj_. 

uration scale in the fundamental representation for both the models can be calculated 
self consistently solving ^^[x,r±^^ = 1/Q'^_p{x,hj_)] = 2(1 — e^^/^). The corresponding 
adjoint saturation scale Q^, relevant for hadronic collisions, is obtained by multiplying 
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Qgp by 9/4. In the range x w 10 ^-10 the behaviour of Qg (see fig.l left) at 6 = 
can be approximated by a function of the form (xq/x)^ with A ^ 0.12 for the b-CGC 
model and A ^ 0.2 for the IP-Sat model. 




Fig. 2. Left: Impact parameter dependence of the probability distribution for inelastic collision. Right: 
Charged-particle multiplicity distributions in |?7| < 1.0 as predicted from saturation models, data points 
are from Ref[18] 

3. Results and discussion 

Multiparticle production in high energy hadronic collisions can be treated self consis- 
tently in the CGC approach. The Glasma flux tube picture [12] predicts [13] that the 
n-particle correlation is generated by the negative binomial distribution {n, k). It 
is characterized by two parameters, the mean multiplicity n and k. At a given impact 
parameter of the collision, the mean multiplicity n « n(b) is obtained by integrating 
eq. 2 over pj^. In the Glasma picture, the parameter k{b) — Ci^c ~ l)Q|5'j_/27r with 
C ^ 0{1) [14]. The quantity Q^S± shown in fig.l (right) is the number of flux tubes 
in the overlap area S± of two hadrons. Convolving P^^{n{b), k{b)) with the probability 
distribution '^^'g'^- for an inelastic collision at 5-fig. 2 (left)-one obtains [6] the n-particle 
inclusive multiplicity distribution as shown in fig. 2 (right). 

Various kinematic variables exhibit scaling with the saturation scale [4,6]. The mid- 
rapidity multiplicity density scales with functional forms like QKs) and Q's{s)/as{Qs) 
whereas a linear functional form seem to provide very good fit to the energy dependence 
of (p±) as shown in fig. 3 [left]. These results are suggestive that is the only scale that 
controls the bulk particle multiplicity. In Ref. [4,7] it has been shown that pj^ spectra in 
p + p collisions exhibit geometric scaling assuming a simple form of Qs- In our case we 
use a scaling variable p±/Qs, where Qs is directly calculated in the IP-Sat model. As 
shown in fig. 3 [right], an approximate scaling below Pa_/Qs < 3 is observed for transverse 
momentum distribution in p + p collision energy ^/s > 540 GeV. Going to lower energies 
we observe systematic deviations from the universal curve. 

In summary, our description of multiplicity distribution successfully describes bulk 
LHC p-l-p data. In particular, we observe that the dominant contribution to multiplicity 
fluctuations is due to the intrinsic fluctuations of gluon produced from multiple Glasma 
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Fig. 3. Left: Scaling of (px) with saturation scale. Right:Scaling of px-spectra plotted in terms ofpx/Qs 
; data points arc from Ref[15-17]. 

flux tubes rather than from the fluctuations in the sizes and distributions of hotspots. 
The Px-spectra in p+p at high energies exhibits universal scaling as a function of p±/Qs- 
The observed scaling indicates that particle production in this regime is dominantly from 
saturated gluonic matter characterized by one universal scale Qs- Ridge like two particle 
correlation structures in Arj — A$ in high multiplicity p+p collisions may provide more 
detailed insight into its properties [19]. 

R.V was supported by the US Department of Energy under DOE Contract No.DE- 
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